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Abstract: 

The cross-section for W + W~ — > W + W~ with arbitrarily polarized W bosons is cal- 
culated within the Electroweak Standard Model including the complete virtual and soft- 
photonic 0(a) corrections. We show the numerical importance of the radiative corrections 
for the dominating polarized cross-sections and for the unpolarized cross-section. The nu- 
merical accuracy of the equivalence theorem is investigated in G(a) by comparing the 
cross-section for purely longitudinal W bosons obtained from the equivalence theorem 
and from the direct calculation. We point out that the instability of the W boson, which 
is inherent in the one-loop corrections, prevents a consistent calculation of radiative cor- 
rections to the scattering of on-real-mass-shell longitudinal W bosons beyond 0(a). 
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1 Introduction 



Gauge-boson scattering provides a direct way to probe both the non-abelian and scalar 
sector of the Electroweak Standard Model (SM). Since in lowest order all gauge-boson 
scattering amplitudes involve only interactions between gauge and scalar bosons, the cor- 
responding cross-sections depend very sensitively on the non-abelian and scalar sector of 
the underlying theory This sensitivity is even enhanced for high-energetic, longitudinally 
polarized, massive gauge bosons owing to the presence of gauge cancellations. For these 
reasons, gauge-boson scattering has found continuous interest in the literature [ [j], ||, H, 3 
since the early years of spontaneously broken gauge theories. The sensitivity of gauge- 
boson scattering to the non-abelian gauge couplings has been investigated, for instance, 
in Ref. [§. 

Gauge-boson scattering can be studied in principle at all high-energy colliders, such as 
pp colliders like the LHC [ ||, e + e~ colliders like the proposed NLC [ 0], or fi + fi~ colliders 
[ ||. These reactions naturally appear as subprocesses, e.g. as final-state interactions in 
hadron or lepton collisions. At high energies (E 3> Mw) the incoming particles radiate 
plenty of gauge bosons. Similar to the well-known Weizsacker-Williams approximation 
for photonic reactions also massive vector-boson scattering at high energies can be ap- 
proximated by convoluting the vector-boson cross-section with the corresponding flux of 
vector bosons. This approximation is known as the equivalent vector-boson method (see 
e.g. Ref. [ |9| and references therein). 

Gauge-boson scattering reactions have been studied at tree level in Ref. [ Q . So far, 
only the leading radiative corrections in the limit of high energies and large Higgs-boson 
masses (E, M H ^> M w ) have been calculated [ |2|, |3|, |J. In most of these calculations the 
scattering of only longitudinal gauge bosons was considered, and the Goldstone-boson 



equivalence theorem [ [K| was used. A complete O(a) calculation was available only for 



the process 77 — > W + W [ 11, 



In a recent paper [ [H| the complete one-loop electroweak radiative corrections to the 
simplest massive gauge-boson scattering process, ZZ — > ZZ, have been discussed. By 
means of this example many characteristic features of massive gauge-boson scattering 
processes could be studied at the one-loop level, such as the treatment of the Higgs- 
boson resonance and the gauge cancellations, which are quantitatively expressed by the 
Goldstone-boson equivalence theorem. 

In this paper we perform a similar study of elastic W + W~ scattering, which is one of 
the dominating gauge-boson scattering processes at all energies. Because the W bosons 
are charged, real bremsstrahlung corrections are required that lead to additional complica- 
tions. Since we are mainly interested in the weak corrections, we include these corrections 
only in the soft-photon limit and omit the radiation of hard photons that depends on the 
experimental situation. Moreover, we restrict our discussion to the case of a relatively 
light Higgs boson so that no Higgs resonance occurs. Such a light Higgs boson is favoured 
by precision experiments [ KM. Furthermore, the treatment of the Higgs resonance re- 



quires to include higher-order corrections. But as we show in this paper, the calculation 
of corrections of order a 2 or higher for the scattering of on-real-mass-shell longitudinal 
W bosons becomes inconsistent owing to the finite width of the W bosons, which enters 
unavoidably via loop diagrams. 
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This paper is organized as follows: After some preliminary remarks about kinemat- 
ics, conventions, and discrete symmetries in Section |2], we discuss the lowest-order cross- 
sections in Section || In Section f| we describe the explicit calculation and the structure of 
the 0(a) corrections. Numerical results are presented in Section |5|. Section |5] is devoted 
to a discussion of problems that occur in the calculation of the radiative corrections to 
W + W~ — > W + W~ beyond O(o). Section |7] contains our conclusions. 

2 Kinematics, notation, and discrete symmetries 

We consider the reaction 

W+(fci, Ax) + W-(A&, A 2 ) -> W+(£;3, A 3 ) + W~(fo, A 4 ) , (1) 

where ki and A, represent the momenta and helicities of the W bosons, respectively. In the 
following, the indices L, T, and U denote longitudinal (Ai = 0), transverse (Ai = ±) and 
unpolarized W bosons, respectively, and definite polarization combinations are labelled 
by a sequence of four letters, e.g. LTLT stands for W^W^ — > W^W^. 

The incoming particles travel along the z axis and are scattered into the x-z plane. In 
the centre-of-mass system (CMS) the momenta and polarization vectors £j(Aj) explicitly 
read 



K = 


(E, 0, 0, -p) , 


— 


(E, —psin9, 0, — pcos9), 


e?(0) = 


(-p, 0, 0, E) /M w , 


eT(0) = 


(p, -Esin6, 0, -Ecos6) /M w 


(±) = 


(0, -1, ±i, 0)/y/2, 




(0, - cos 9, =Fi, sin^) /V2, 




{E, 0, 0, p) , 


= 


(E, psin9, 0, pcos9) , 


4(0) = 


(-p, 0, 0, —E) /M w , 


#•(0) = 


(p, Esin9, 0, Ecos9) /M w , 


= 


(0, 1, ±i, 0)/y/2, 


^*(±) = 


(0, cos 9, =Fi, - sin 9) /V2 



(2) 



in terms of the energy E of the W bosons, their momentum p = a/ E 2 — M w , and the 
scattering angle 9. The Mandelstam variables are defined as 

s = (h + A; 2 ) 2 = 4E 2 , 

t = {kx - hf = -Ap 2 sin 2 9/2, (3) 
u=(k 1 - k A ) 2 = -Ap 2 cos 2 9/2 . 

The polarized differential cross-section is obtained from the invariant matrix element 
M as 

dcr\ 1 , . , l2 



- 



dQ J x x , , 64tt 2 s 
/ A1A2A3A4 



|-M AlA2 A3A 4 | , (4) 



and the unpolarized cross-section results from averaging over the polarizations of the 
incoming and summing over the polarizations of the outgoing particles, 

U1J / unpol y A 1 ,A 2 =-1 A 3 ,A 4 =-1 V J A!A 2 A 3 A 4 
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More generally, the correct average is obtained by multiplying with 1/3 for each unpolar- 
ized W boson and by 1/2 for each transverse W boson in the initial state. 
The integrated cross-section is defined by 

/•360° rl80°-e cut j 

a = / dtp d6 sm6 — , (6) 

dll 



cut 



where 8 cut denotes an angular cut which is set to 10° in the numerical calculations. 
Discrete symmetries relate different helicity amplitudes. CPT symmetry implies 

Mx^iXsX^ = M\ 3 \ 4 x 1 x 2 . (7) 

Since we neglect quark mixing, also CP is conserved and as a consequence, the helicity 
amplitudes are related as follows 

MxiXiXa^ — M-X 2 -Xi-Xi-X 3 ■ (8) 

These symmetries reduce the number of independent polarized amplitudes from 81 to 27. 

3 Lowest-order cross-sections 

3.1 General properties 

In lowest order, seven Feynman diagrams exist for W + W~ — > W + W~ as shown in 
Fig. [1]. Owing to charge conservation there are no w-channel diagrams, and therefore the 
cross-section is expected to have a forward-backward asymmetry. The resulting tree-level 
amplitude reads 

M BoTn = -e 2 (± + ^^—) Ml -e 2 (l+&-L-)M l 



2 2 M 2 i 2 M 2 i W 

s 2 w s 2 w s-M| s 2 w t-Ml 
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The matrix elements M.1-M.5 correspond to the diagrams #l—#5 in Fig. [I], respectively, 
and contain the complete kinematical information and polarization dependence, 



-Mi = (£1 • £2) (£3- el)(t-u) 

+ 4{( £l • e 2 ) [{el ■ fa) (e* 4 ■ fa) - (e* 3 ■ fa) (e* 4 ■ fa)] 

+ [el ■ el) [(£1 • fa){e 2 ■ fa) - ( £l -fa){e 2 -fa)] (10) 

- (e 2 • fci) [(£1 • £3) (£4 • fa) - (£1 • £4) (£3 • fa)] 

- (£1 • k 2 ) [(e 2 ■ el) (el ■ fa) - (e 2 • £3) (£4 • fa)] } , 
M 2 = {e 1 -el){e 2 -el){s-u) 

- 4{( £l • e* 3 ) [(e 2 • fc 3 )(eS • fc 2 ) + (£2 ■ fc 4 )(e2 • fc 3 )] 

+ (£2 • el) [( £l ■ A; 2 )(e* ■ fe 4 ) - ( £l ■ fc 4 )( £ * ■ fc 2 )] (11) 

- {e x ■ fa) [{el ■ el)(e 2 ■ fa) + (e 2 ■ e*)(ej • fc 2 )] 

- (£3 ■ fa) [(£1 • £ 4)(£ 2 • fa) + (£1 • e 2 )(e\ ■ fa)] J , 

M 3 = (£! • £ 2 )(£* • £ 4 ) + ( £l • £^(£ 2 • ej) - 2{e x ■ e* 4 )(e 2 ■ e* 3 ) , (12) 
M A = (e 1 -e 2 ){el-el), (13) 

M 5 = (£!• £3)(£2"£l). (14) 

For equal polarizations in the initial and final state, i.e. for Ai = A3 and A 2 = A 4 , the 
lowest-order amplitude diverges for small momentum transfer t owing to the t-channel 
photon-exchange diagram. In fact, in this limit the lowest-order cross-section turns into 
the Rutherford cross-section oc 1/t 2 oc l/[psin(6'/2)] 4 . If either Ai = A 3 or A 2 = A 4 , the 
cross-section has at most a 1/t singularity for small t, and for Ai 7^ A3 and A 2 7^ A 4 it is 
regular. 

Each longitudinal W boson enhances the individual matrix elements Aii-Ai?, by a 
factor -y/s/Mw. The cancellations guaranteed by gauge symmetry prevent the cross- 
section for longitudinal WW scattering from grossly violating unitarity at high energies. 
The gauge cancellations are demonstrated graphically in Fig. |2] for the total cross-section. 
The cross-section calculated from individual diagrams involving three- (a) or four-gauge- 
boson couplings (b) grows with s 3 for large s. Adding up all pure gauge-boson diagrams 
the leading terms cancel resulting in a cross-section oc s. Including the Higgs-exchange 
diagrams, further cancellations take place such that the SM cross-section drops as 1/s. 
The cancellations are already quite substantial at a few hundred GeV. 

The unitarity cancellations ensure that all polarized cross-sections decrease at least as 
1/s at high energies. The cross-sections with an odd number of longitudinal W bosons 
are suppressed by an additional factor M^/s. The behaviour of various polarized cross- 
sections at high energies is listed in Table [IJ. 
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Polarization 


o"Born = vra 2 / ' • (entry) 


LLLL 


1 ( (M 2 + M|) (2Mf L c + M 2 c) c (75 - 26c 2 - c 4 ) \ 




4 4M^ 1 484(1-^) J 


LLLT 
= LLTL 
= 2 LTLL 
= 2 TLLL 


1 ( (M 2 + 4M 2 , + M 2 ) 2 L C 
s 2 I 4M 2 , 

c[3 (M 2 + 2(3M 2 , + M 2 )) 2 - (6 - c 2 )(2M^ + M 2 ) 2 ] ] 




12M 2 , J 


LLTT 
= 4TTLL 


lc(3 + c 2 ) 


LTLT 
= TLTL 




LTTL 
= TLLT 


1{ [M 4 + M| + 2M 2 (8M^ + M|)]L C + ° (5 |^f^ 
| c (8M 2 , + Ml) [18M 2 + (3 + c 2 )(8M^ + M 2 )] ] 


) 




LTTT 
= TLTT 
= 2 TTLT 
= 2 TTTL 


-^M 2 ,(29c + 3c 3 + 20L c ) 


TTTT 





Table 1: Leading terms for polarized integrated lowest-order cross-sections in the high-energy 
limit. Here c = cos# cu t, L c = ln^xf, and the relations between different polarizations in the 
left column are to be read, for instance, as olttt = 2<tttlt> where the numerical prefactors 
originate from the different spin averages. 



3.2 Numerical results 

All following numerical results are obtained for an angular cut-off of 8 cnt = 10° and for 
a Higgs-boson mass Mh = 100 GeV. This value is within the range allowed by precision 



data, which favour a light Higgs boson [p!4[. Since it lies below the threshold for W + W _ — > 
W + W~ no Higgs resonance occurs. 

In Fig. ^| the integrated lowest-order cross-section is displayed for various W polariza- 
tions. The curves nicely reflect the high-energy behaviour of the analytical expressions in 
Table |l|. The cross-sections for the polarizations LLLL, TTTT, LTLT, and LTTL behave 
as 1/s at high energies, all others are suppressed by additional factors M^/s. Owing 
to the 1/t pole, the cross-sections with equal initial- and final-state polarizations, LLLL, 
TTTT, and LTLT, are enhanced by a factor 1/(1 — cos 2 9 cut ) ~ 33 for Q cut = 10° and are 
therefore dominant. For these polarizations the contribution of the backward hemisphere 
to the integrated cross-section is less than 1% for energies above 1 TeV and at most 3% for 
lower energies. Note that the size of the corresponding integrated cross-sections strongly 
depends on the angular cut, cut , because of the t-channel singularity. 
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Figure 2: The cross-sections for longitudinal gauge-boson scattering resulting from subsets of 
the tree-level diagrams: (a) diagrams involving only three-gauge-boson couplings, (b) diagram 
involving only four-gauge-boson couplings, (c) diagrams involving Higgs bosons. 



or 



orn 10^ 
10 2 
10° 
lO" 2 h 
10 
10 
10 



1 1 1 


1 1 










- v " -O" - 








M H = 100 GeV 
l l l 


1 1 



TTTT 
UUUU 
LTLT 
LLLL 
LLTT 
LTTT 
LLLT 



200 500 1000 2000 5000 10000 y^/GeV 

Figure 3: The integrated lowest-order cross-sections for various polarizations. 
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Figure 4: The differential lowest-order cross-sections for various polarizations at different ener- 
gies. 



In Fig. [| the differential lowest-order cross-sections are shown for a/s = 200 GeV, 
1 TeV, and 5 TeV for the dominating polarizations LLLL, LTLT, LLTT, TTTT, and for 
the unpolarized case UUUU. For all these polarizations except LLTT the curves display 
the strong forward peak resulting from the t-channel pole. The figures show that the 
main contributions to the unpolarized cross-section in backward direction come from 
completely transversally polarized bosons. The sharp dips in the longitudinal cross-section 
correspond to zeros in the amplitude. To be precise, the longitudinal cross-section is a 
rational function in t and has four such zeros. For yfs < 274 GeV two of these zeros and 
for -^/i > 274 GeV one of them lie in the physical region. 
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4 O(et) corrections 



4.1 Calculational framework 

We have performed the calculation of the 0(ct) radiative corrections in 't Hooft- 
Feynman gauge both in the conventional formalism and in the background-field formalism. 
Ultraviolet (UV) divergences are treated within dimensional regularization. We use the 



on-shell renormalization scheme [ [H| , following the formulation worked out in Ref . [ 16] 



for the conventional formalism and in Ref. [ O for the background-field formalism. In 



the conventional formalism the field renormalization is fixed such that no external wave- 
function renormalization is needed. In the renormalization scheme introduced in Ref. f IT? 



for the background-field method the field renormalization is determined by gauge invari- 
ance, and a non-trivial external wave-function renormalization is required, as explicitly 
described in Ref. [ [18 . 



Our calculation is based on the second of the two methods described in Ref. [ |13| 



The Feynman diagrams are generated with FeynArts [ [L9|] . The resulting amplitudes are 
algebraically simplified with a combination of Form [ |2(| and Mathematica and are then 



converted into a Fortran program. The tensor integrals are numerically reduced to scalar 



integrals, which are evaluated using the FF package [ [21 
4.2 Inventory 

At one- loop level, roughly 1000 diagrams contribute to W + W~ — > W + W~. The 
diagrams can be classified into self-energy corrections, vertex corrections, and box correc- 
tions. The self-energy corrections and the vertex corrections can be further divided into 
s- and t-channel contributions, the box corrections involve in addition u-channel contribu- 
tions. The t-channel contributions can be obtained from the s-channel contributions via 
crossing, i.e. via the interchange of the outgoing W~ boson with the incoming W + boson. 

In the following we list the Feynman graphs in the conventional formalism in 't Hooft- 
Feynman gauge. To avoid presenting hundreds of diagrams, similar diagrams have been 
combined. Each combination of particle labels on the internal lines corresponds to an extra 
diagram, where those have to be omitted that violate charge or lepton-number conserva- 
tion or involve one of the following vertices: yyHH, 77XX, 'yZHH, ^ZxXi W + W~\H , 
xHH, Z\X) ZHH, 7XX> 1XH, ^HH, juu, which do not exist in the SM. For example, 
diagram #1 of Fig. |5| represents a 0-loop contribution to the photon self-energy and to 
the 7Z and Z7 mixing, and 0-, x~i an d H-loop contributions to the Z-boson self-energy. 
If the charge flow in the internal lines is not indicated, it can take both directions (e.g. 
in diagram 7^5 in Fig. [5]) or it is determined from the charges of the external W bosons 
(e.g. in diagram #2 in Fig. |6]). In addition, fermionic diagrams are shown only for the 
first fermion generation. 

The s-channel self-energy corrections are shown in Fig. EL They consist of insertions 
of the Z-boson and photon self-energies and the 7-Z-boson mixing energy into diagram 
#1 of Fig. |l] and Higgs-boson self-energy insertions into diagram $4 of Fig. [I]. 

The s-channel vertex corrections consist of corrections to the 'yWW, ZWW and 
HWW vertices in diagrams #1 and $4 of Fig. |l[ The corrections to the final-state 
vertices in these diagrams are displayed in Fig. 0. 
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The box diagrams are shown in Figs. [7|, || and [3]. While the t-channel boxes are 
obtained from the s-channel boxes via crossing, the w-channel boxes are different owing 
to the different charge flow. In particular, there are no fermionic w-channel diagrams. For 



this reason, the Landau singularity that appears in the corrections to ZZ — > ZZ [ [L3[ does 
not show up in W + W~ — > W + W~. 

Apart from the virtual corrections, one has to take into account real soft-photon 
emission from the external W bosons in order to cancel the infrared (IR) divergences in 
the charge form factor. The IR divergences are regularized by a infinitesimal photon mass 
A. In the soft-photon limit the cross-section for real photon emission is given by 

(£) =(£) 4* M 



V dfi / so ft Vdfi/Born 
with the soft-photon correction factor 



'soft 



e 2 f d 3 k A ±Q i Q j {h ■ k 



E 



(2t) s J^ae 2*d ^ (k ■ k)(kj ■ k) 



fc = V /fe2 +A 2 



(16) 



Here AE is the maximum energy of the emitted photons, Qi = ±1 is the charge quantum 
number of the ith W boson, and the sign in front of the product QiQj is '+' if particles i 
and j are both either incoming or outgoing and '— ' otherwise. 

The basic integrals needed for the soft-photon factor have been worked out, for in- 
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#29 #30 #31 #32 

Figure 6: The s-channel vertex diagrams for the final-state vertex. 
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Figure 7: Genuine box diagrams that appear in the s, t, and u channel. Each diagram in (a) 
represents actually four diagrams with different charge flows shown in (b), i.e. one s-channel, 
one t-channel, and two u-channel diagrams. 



stance, in Ref . [ |22] . For the case of W + W — > W + W the soft-photon factor reads 

2a f 2AE E E-p 

<Wt = S 2 In — 1 In — — — 

7T [ A p E + p 



-{s- 2Ml) X(s + y/s(s - 4M^)) 

- (2M^ - t) j(4M^ - i + v/t(t-4M^)) 

+ (2M^ - it) j(4M^ - u + > /«(u-4Af^))J ! 



(17) 
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\ X,B,<p w 



Figure 8: Further s-channel box diagrams. For every diagram there exists a t-channel analogue 
(obtained by exchanging the lower left and upper right leg). 





Figure 9: Further -u-channel box diagrams, 
where the function X(x] is denned as 
X{x) 



2{x-2M^) / 21n 2AE In x-2M^ 



x(x - 4M^ 



+ Li 2 1 - 



A 



2M 2 
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x 



(18) 



4.3 Leading corrections 

At low energies (E ~ M w ) the electroweak radiative corrections are dominated by the 
universal corrections associated with the running of a and the p parameter. These effects 
are typically at the level of 10%. In this energy region no large logarithmic corrections 
associated with collinear photons appear in W + W~ — ► W + W~ because no light external 
particles are involved. 
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At high energies (E ^> Mw) and for a light Higgs boson, the radiative corrections 
are dominated by corrections of the form (at/ir) \n(\x\/M^) ln(|y|/M^), x,y — s,t,u, 
which originate from approximate IR and collinear singularities in the vertex and box di- 
agrams.Q These terms are are of the order of 10 and 20% at 1 and 10 TeV, respectively. At 
even higher energies these logarithms render the one-loop results useless and perturbation 
theory breaks down. In order to obtain sensible predictions for exclusive cross-sections 
beyond 10 TeV some kind of resummation of these logarithms is required. When consider- 
ing more inclusive quantities and taking into account W-boson, Z-boson, and Higgs-boson 
bremsstrahlung, the (a/vr) ln(|a;|/M^) \n(\y\/M^) terms should drop out leaving leading 
terms of the form (cx/tt) ln(|x|/M^). 

If the Higgs boson is very heavy, additional large corrections appear. These correc- 
tions are, in particular, important for the scattering of longitudinal W bosons and have 
been extensively discussed in the literature [ 0, ^3). In the limit M| ^> s ^> M^, 



terms of the form (a/nsfa)(\x\/M$ f )hi(M&/\y\), x,y = s, t, u, which have been cal- 
culated in Refs. [ |2|, |24]|, dominate the relative corrections. For s ^> ^> M^, 
on the other hand, the leading relative correction is provided by terms of the form 
(a/7rs^ v )(M^/M^ v )ln(\x\/M^ L ). In this limit also the two-loop corrections are available [ 
p5| . The complete one-loop corrections proportional to M^/M^ to longitudinal W-boson 
scattering for s, M\ 3> can be found in Ref. [ [J. Finally, the leading logarithmic 
corrections of the form (a/ns^) W-^h/^w) ^ n ^ ne nm ^ -^h ^ s > f° r the scattering 
of arbitrarily polarized W bosons are given in Ref. [ ^] . 

4.4 Checks of the calculation 

Several cross-checks were performed to ascertain the correctness of our calculation. In 
particular, we checked 

• the algebraic simplification of the Feynman diagrams by comparing the Form results 
with those obtained from FeynCalc [ ^7]] for selected diagrams (analytically), 

• the numerical evaluation of the Feynman diagrams by comparing the evaluation of 
selected diagrams in Fortran and in Mathematica (numerically), 

• the evaluation of the scalar and tensor one-loop integrals by comparing the results 
of the FF routines with own routines (numerically), 

• the gauge independence of the results by comparing the results in conventional 
formalism and background- field formalism (numerically), 

• the UV finiteness by proving independence of the scale parameter \x of dimensional 
regularization (analytically), 

• the IR finiteness by testing independence of the regulator mass A (numerically). 

Numerical computations were done using extended precision (quadruple precision) with 
a mantissa of approximately 33 digits. The relative deviation between numerical results 

1 For the process e + e _ — * W + W~ these leading terms can be read off from the high-energy expansion 
given in Ref. [ E3| . 
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from different calculational methods was less than 10~ 16 . A double precision calculation 
is in general sufficient up to about 10 TeV. For higher energies the gauge cancellations 
become so large that double precision does not give reliable results any more. 

5 Numerical results for the O(ct) corrections 

5.1 Input parameters and definition of the corrected cross-section 



For the calculations the following parameter set is used [ [28 1: 

a' 1 = 137.0359895, M z = 91.188 GeV, M w = 80.401 GeV, 

m e = 0.51099906 MeV, m u = 47.0 MeV, m d = 47.0 MeV, 
m M = 105.658389 MeV, m c = 1.55 GeV, m s = 150 MeV, 
m T = 1771.1 MeV, m t = 175.5 GeV, m b = 4.5 GeV. 



(19) 



The masses of the light quarks are adjusted such that the experimentally measured 
hadronic vacuum polarization is reproduced [ |29j. As before, a Higgs-boson mass of 
Mh = 100 GeV is used, and the integrated cross-sections are evaluated for an angular cut 

0cut = 10°. 

All cross-sections discussed in the following are calculated by replacing the matrix 
element squared in (|J) by 

\M\ 2 - l-MBoml 2 (1 + 4oft) + 2 Re(M* Born 5 M) , (20) 

where 5Ai is the sum of the one-loop Feynman diagrams, and 5 so f t the soft-photon cor- 
rection factor (|T6|) . Note that this yields the cross-section including all 0(a) corrections 
but neglecting all 0(a 2 ) corrections. 

In the soft-photon approximation the corrections depend on the soft-photon cut-off 
energy AE via logarithms of the form \n(AE/E). For small AE, where the soft-photon 
approximation is valid, these terms give rise to large corrections, which are, however, 
cancelled if hard-photon bremsstrahlung is taken into account. In order to avoid artificially 
large soft-photon effects, we have chosen to discard all terms involving \a(AE/E) in the 
corrections.^ The elimination of the cut-off-dependent terms can be simply achieved by 
choosing AE = E. The corrections defined in this way can be viewed as a suitable 
measure of the weak corrections, which cannot be separated from the electromagnetic 
corrections on the basis of Feynman diagrams in a gauge-invariant way [ 1 1 1 . 



5.2 Total cross-sections 

In Fig. [TO, the integrated cross-sections in lowest order and including 0(a) corrections 



as well as the relative O(a) corrections are plotted for the dominant polarizations LLLL, 
LTLT, LLTT, TTTT, and UUUU. For small energies the corrections are positive and 
roughly 10%. With increasing energy the corrections become large and negative and 
reach the order of the lowest-order cross-sections around 10 TeV. The large corrections at 
high energies originate from logarithms of the form a/7rln (s/M^), which are additionally 
enhanced by numerical factors. 

2 Of course, these terms have to be reinstalled in a complete calculation involving hard-photon 
bremsstrahlung. 
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Figure 10: Upper plot: The integrated cross-sections for the dominant polarizations in lowest 
order (Born) and including the 0(a) corrections (corr). Lower plot: Relative corrections for the 
same polarizations. 



5.3 Differential cross-sections 

The lowest-order and corrected differential cross-sections for the dominant polariza- 
tions LLLL, LTLT, LLTT, TTTT, and UUUU are shown in Fig. O for ^/s = 200 GeV, 



1 TeV, and 5 TeV. The corrections are small at low energies for all scattering angles but 
get large at high energies, in particular in the backward direction. This is due to the fact 
that the contributions of the w-channel box diagrams to the cross-sections involve terms 
that behave as l/\u\ for <C \u\ <C s and approach a constant of the order of 1/M^, 
for \u\ <C Myj, while the lowest-order cross-sections are of order 1/s for \u\ <C s. Hence 
the relative corrections in the backward direction are typically enhanced by a factor s/\u\ 
compared with the corrections in forward direction or to the total cross-section and can 
easily reach 100% at several TeV. Since |<5.M| 2 is missing in (^0|), the cross-section can 
formally become even negative in this case. By including |<5.M| 2 , the quality of the predic- 
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Figure 11: The differential cross-sections for the dominant polarizations at y/s = 200 GeV, 
1 TeV, and 5TeV in lowest order (Born) and including the 0(a) corrections (corr). Note that 
in some regions of phase space the corrected cross-sections become formally negative and are 
therefore not visible in these logarithmic plots. 

tions can in general be improved, but for W + W~ — > W + W~ this leads to problems, which 
are discussed in Section ||. Here, we ignore |<5.M| 2 and the resulting wrong predictions 
for high energies and large scattering angles, since the cross- sect ions are small in this 
region, and, as has been stated above, the total cross-section is dominated by the forward 
hemisphere. The large corrections in the backward direction do not indicate a breakdown 
of perturbation theory because no additional enhancement factors s/u should appear in 
higher orders. 
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5.4 Accuracy of the equivalence-theorem predictions 



The calculations of enhanced corrections to longitudinal gauge-boson scattering [ ^ 
usually have been done with the help of the equivalence theorem (ET) [[RJ. When using 
the ET in higher-order calculations one has to be careful to include all correction factors 
that result from renormalization and amputation 
the BFM 



30 . This is particularly easy within 



IS 



where the matrix elements for external longitudinal vector bosons are 
directly obtained from the amputated Green functions with the corresponding would-be 
Goldstone-boson fields multiplied with the wave-function renormalization constants of the 
gauge bosons. 

We calculated the matrix element for tp + tp~ — > (f +i p~ to one-loop order in the BFM 
and from this the ET prediction for the cross-section for Wj^W^ — > Wj^W^ 



da\ 



1 



QAtt 2 s 



I \/t'P + 'P -^^v 



Born 



[i + s, 



soft J 



+ 2 Re 



{(M 



Born 



TV 



)*5M 



[21) 



Since y? + y2~ — > tp + <{>~ and (p + ip~ — > Lp + (p~-y are no physical processes, the cross-section 
is not well-defined at finite energies. The ET is valid only for high energies and 
guarantees ( |21| ) to provide a sensible cross-section only in the high-energy limit, i.e. up 
to terms of order Myy/E. In fact, at finite energies fl2l|) is gauge- dependent and the IR 



singularities between the virtual and real corrections do not cancel exactly. In the analytic 
results, which are obtained for an infinitesimal photon mass A, terms proportional to 
(M-w/E) ln(A 2 /s) and gauge-dependent terms proportional to (M w /E) remain. In order 
to avoid an enhancement of the IR-singular terms, we have set the photon mass equal 
to the energy of the W bosons, A = E, in the numerical evaluation. This ensures that 
the artificial IR-singular logarithms do not lead to enhanced contributions. Since also 
the soft-photon correction factor corresponding to (p + tp~ — > y? + y?~ is gauge- dependent we 
have chosen to use directly the soft-photon correction factor for W + W _ — ► W + W~. We 
stress that despite the gauge-dependences and the IR-singularities contained in (|2l|) , our 
results are valid at the same level as all ET results, i.e. up to terms of order M-w/E. 

In order to improve the accuracy of the ET calculation we combine the lowest-order 
cross-section for Wj^W^ — > Wj^W^ with the 0{a) corrections from <fi + tp~ — > y? + y?~ as 
follows 



da\ 

mixed 



64vr 2 s 



i » iW+wr- 

•^Born 



♦W+Wl|2 



K 'Born 



soft 



+ 2Re<! {M v+Lp 



Born 



)*5M 



(22) 



In this way the error originating from the use of the ET enters only via the 0(a) correc- 
tions. 

The accuracy of the ET predictions for the scattering of longitudinal W bosons is inves- 
tigated in Fig. [12] and Table 0. We show the lowest-order cross-section for <p + (f~ — > 
normalized to the lowest-order cross-section for W^Wr — > W^Wr and the corrected 
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Figure 12: Ratio of the ET predictions for W^Wj^ — > WjW[ and the direct calculation. "Born" 
and "corr" label the ratios of the tree- level and 0(a)-corrected cross-sections, respectively. For 
the curve marked "mixed" the ET is used only for the corrections [cf. (f^)]. 





Born 


o — > tp + (p' 
Tl -> W+\ 
corr 


mixed 


500 GeV 


-6.4% 


-2.7% 


3.2% 


ITeV 


-3.0% 


0.8% 


3.7% 


2TeV 


-1.0% 


0.3% 


1.3% 


5TeV 


-0.2% 


-0.6% 


-0.4% 


lOTeV 


-0.0% 


-0.2% 


-0.1% 



Table 2: Relative deviations of the ET predictions for Wl W L — ► W L . 



cross-section calculated from the ET according to and according to the improved 
formula (0), both normalized to the corrected cross-section for Wj^W^ — > W^W^- For 
y^s > 500 GeV the deviations are all below 7% and become smaller at higher energies. 
Above about 2TeV they are less than 1%. At high energies, where the ET is valid, the 
approximation (f22|) is in fact better than the pure ET prediction. 

6 Difficulties beyond O(ct) accuracy 

In Section |5] we found that the 0(a) corrections become of the order of the lowest- 
order cross-section for high energies. In this case higher-order corrections are required in 
order to arrive at a meaningful prediction. If the importance of the 0(a) corrections is 
due to a suppression of the lowest-order cross- sect ion, as it happens for W + W~ — > W + W~ 
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at large scattering angles and high energies, it is appropriate to include the square of the 
0(a) corrections, i.e. to calculate the cross-section via 

M = ^ M ^ + SMf - (23 > 

Because the interferences of the genuine 0(a 2 ) matrix element with the lowest-order 
matrix element are suppressed with the lowest-order cross-section, this gives a reasonable 
(leading-order) prediction, which is usually sufficient as the cross-section is suppressed. 
This recipe has for example been applied in the case of ZZ — ► ZZ [ p~3f] . 

Since for W + W~ — > W + W~ this approach leads to severe problems, it was not used 
in this paper. 

The first problem is due to the charge of the W boson. It gives rise to photonic 
corrections that render the matrix element for the process W + W~ — > W + W~ IR-singular. 
In order to cancel the IR singularities one has to add real soft-photon radiation. The 
cancellation of the IR singularities takes place order by order. Including \8M. \ would thus 
require to take into account all virtual and real 0(a 2 ) corrections. Since the IR-singular 
corrections are proportional to the lowest-order matrix element and thus suppressed, an 
approximate result could nevertheless be obtained by including the real 0{a) corrections 
appropriately when squaring the matrix element. The error introduced by this effective 
treatment of the real corrections is of the order of the neglected 0(a 2 ) corrections. 

The second, more severe difficulty is closely related to the instability of the W boson. 
When calculating the one-loop diagrams for WlW l — > W^W L one finds that not all 
terms proportional to s cancel. The left-over terms proportional to s in the unrenormalized 
amplitude can be expressed in terms of the transverse W self-energy: 

;jU« _ ^ cos 2 g S ^T( M w) m) 
d ^O(s) ~ 2 COS 2M 2 M 2 • W 
6 W z 1V1 W 1V1 W 

The counter terms, in fact only the W-mass counter term 5M^ = Re£^(M{y) (and 
the derived counter term for the weak mixing angle Sc 2 ^), yield a similar contribution 
with opposite sign. When adding these contributions, the real part of S-M^™' is exactly 
cancelled but the imaginary part survives. As a consequence the imaginary part of the 
one-loop matrix element for Wj^W^ — > Wj^Wj^ grows with s and therefore would violate 
unitarity if taken into account. Because the lowest-order matrix element is real as long as 
no Higgs resonance is present, this does not affect the 0(a) corrections to the cross-section. 
However, in higher orders or if the Higgs-boson width has to be included the imaginary 
part becomes relevant and the perturbative evaluation of corrections to W^W^ — > W^W^ 
for on-real-mass-shell external W bosons breaks down. This statement holds for the usual 
approach where the external W bosons are treated as stable particles with a real mass, 
i.e. where the corresponding momenta squared are set equal to the real W-boson mass 
squared. 

The reason for the breakdown becomes apparent by noting that the unitarity- violating 
term can be written as 



A-koi 2 9 s iImSJ l 5 / (iW w ) .4vra , 9 s i w 

„ _^ = 1 „ _ , (25) 
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i.e. it is proportional to the (lowest-order) decay width of the W boson. This result 
shows that the problem is due to the instability of the W bosons, which is inherent in the 
one-loop calculation. It points to the fact that S-matrix elements for unstable external 
particles are not well-defined. 

From the previous considerations it is clear that the unitarity-violating terms would 
cancel if the W-boson-mass counter term were complex, i.e. if renormalization would also 
compensate the imaginary part of the W self-energy. Since the bare W-boson mass is 
real, this can only be achieved if the renormalized W-boson mass (and the renormalized 
weak mixing angle) become complex. 

Let us consider the matrix element in this renormalization scheme including 0(a) 
corrections. In the one-loop part and in the counter-term part, which are of order a, the 
renormalized parameters can be replaced by bare parameters and all unitarity-violating 
terms cancel out. In the tree-level matrix element, however, the complex renormalized 
parameters lead to unitarity-violating terms (if it is calculated for kf = M^). This reflects 
the fact that a change in renormalization scheme in an 0(a) calculation cannot eliminate 
O(a) corrections but only reshuffle them. However, if we evaluate the matrix element for 
complex external momenta squared, kf = — iM-wIV, all unitarity-violating terms 
cancel in the tree-level part, and the loop and counter-term parts are not changed (up to 
terms of order a 2 ) with respect to the evaluation with real external momenta squared, 
kf = M^j. As a result, we find that the matrix element for WJWl — > WJWl including 
O(a) corrections respects unitarity if the squares of the external momenta are equal to 
the physical complex masses, but violates unitarity for real "on-shell" external momenta. 
This result is independent of the renormalization scheme, and holds, in particular, in the 
usual on-shell scheme. Our findings are in accordance with the fact that gauge invariance 
and unitarity of S'-matrix elements require the momenta of the external particles to be on 
the mass shell, i.e. on the complex mass shell for unstable particles. They indicate that 
the natural generalization of S'-matrix elements for unstable particles are the multiple 
residues at the physical (complex) poles of the external particles (compare Ref . [ ^TJ ) . 

The approach described above requires to consider matrix elements for complex exter- 
nal momenta squared. However, the momenta of physical particles are always real, even 
if these particles appear only as resonances in some process. Since it is not clear how 
to relate the physical real momenta to appropriate complex momenta we have not pur- 
sued this approach any further. The definition of S-matrix elements for unstable external 
particles remains an unsolved problem. 

When calculating the corrections to W^W^ — > W^W^ via the ET, no unitarity- 
violating terms appear. The mismatch between the direct calculation and the ET calcu- 
lation, i.e. the apparent violation of the ET, can be understood by looking at the derivation 



of the ET [ 18 . The ET results from Ward identities for connected Green functions by 
amputation and putting external fields on their mass shell. Thus, the validity of the ET 
demands external W bosons to be on their mass shell, i.e. on the complex mass shell for 
unstable particles. For real external on-shell momenta the ET is violated beyond 0(a) 
by terms proportional to the W-boson decay width that are enhanced by factors E/Mw 
originating from longitudinal polarization vectors. The failure of the ET is just another 
symptom of the lack of a proper description of S-matrix elements for unstable particles. 
The only fully consistent way to calculate corrections to longitudinal W-boson scatter- 
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ing beyond 0(a) seems to be the evaluation of a physical process, e.g. e + e~ — > W + W~ — > 
W + W~ — > 4 fermions, or at least the resonant contributions to this process. This is 
beyond the scope of this paper. 

7 Conclusions 

The scattering of the weak gauge bosons is an extremely useful tool for the investiga- 
tion of the non-abelian gauge sector and the symmetry-breaking sector of the electroweak 
interaction. To lowest order these reactions involve only interactions between gauge and 
scalar bosons and, therefore, depend very sensitively on these couplings. Owing to its large 
cross-section, a particularly interesting process in this class is the scattering of charged 
W bosons, W+W- -> W + W~. 

This process is known in the literature including only the leading corrections in the 
limit of high energies and large Higgs-boson masses. In this paper we have extended 
these results by calculating the complete virtual and soft-photonic 0(a) corrections to 
the amplitudes for W + W~ — > W + W~ as predicted by the electroweak Standard Model. 

The cross-sections with equal initial- and final-state polarizations exhibit a 1/t 2 singu- 
larity resulting from the t-channel photon-exchange diagram and are therefore dominant. 
The corresponding 0(a) corrections are typically of the order of 10% for energies below 
roughly 2 TeV, grow with energy, and reach the order of the lowest-order cross-section at 
several TeV. The corrections become particularly large in the backward direction, where 
the cross-sections are relatively small. At energies beyond 10 TeV the one-loop calculation 
becomes unreliable and a resummation of the leading corrections or the inclusion of real 
massive boson emission is required. 

The predictions of the equivalence theorem, which relates the scattering amplitudes 
of longitudinally polarized gauge bosons with those of the unphysical Goldstone-boson 
fields, have been compared with the direct calculation at the level of the 0(a) corrections 
and for a light Higgs boson. For energies above 500 GeV the deviations are below 7%, 
above about 2 TeV they are less than 1%. 

We have found that the calculation of radiative corrections to the scattering of on-real- 
mass-shell longitudinal W bosons, i.e. to W^Wl — > WJWl in the usual sense, becomes 
inconsistent beyond 0(a). More precisely, the imaginary part of the one- loop matrix 
element gets a unitarity-violating contribution proportional to the W-boson width. This 
contribution does not affect the 0(a) corrections to the cross-section as long as no Higgs 
resonance appears. It becomes problematic as soon as the imaginary part of the matrix 
element enters an observable, i.e. if 0(a 2 ) corrections to the cross-section or the finite 
width of the Higgs boson are included. This result points to the fact that S'-matrix 
elements for unstable external particles are not well-defined. For a consistent evaluation, 
the instability of the W bosons has to be taken into account properly. This seems to be 
only possible when considering a full process with stable initial- and final-state particles 
that involve W^W^ — > WlW[ as an enhanced contribution. 
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